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Abstract. Since the expected rates for neutralino-nucleus scattering are expected to 
be small, one should exploit all the characteristic signatures of this reaction. Such 
are: (i) In the standard recoil measurements the modulation of the event rate due to 
the Earth's motion, (ii) In directional recoil experiments the correlation of the event 
rate with the sun's motion. One now has both modulation, which is much larger and 
depends not only on time, but on the direction of observation as well, and a large 
forward-backward asymmetry, (iii) In non recoil experiments gamma rays following 
the decay of excited states populated during the Nucleus-LSP collision. Branching 
ratios of about 6 percent are possible (iv) novel experiments in which one observes the 
electrons produced during the collision of the LSP with the nucleus. Branching ratios 
of about 10 per cent are possible. 

1 Introduction 

It appears without doubts that dark matter constitutes about 30 % of the energy 
matter in the universe. The evidence comes from the cosmological observations 
PP, which when combined lead to: 

^2b = 0.05, ncDM = 0.30, Ha = 0.65 

and the rotational curves pi. It is only the direct detection of dark matter, 
which will unravel the nature of the constituents of dark matter. In fact one 
such experiment, the DAMA, has claimed the observation of such signals, which 
with better statistics has subsequently been interpreted as modulation signals [3]. 
These data, however, if they are due to the coherent process, are not consistent 
with other recent experiments, see e.g. EDELWEISS and CDMS 
Supersymmetry naturally provides candidates for the dark matter constituents. 
In the most favored scenario of supersymmetry the LSP can be simply described 
as a Majorana fermion (LSP or neutralino), a linear combination of the neutral 
components of the gauginos and higgsinos 0-|Hl- We are not going to address 
issues related to SUSY in this paper, since thy have already been addressed 
by other contributors to these proceeding. Most models predict nuclcon cross 
sections much smaller the the present experimental limit as < lO^^pb for the 
coherent process. As we shall see below the constraint on the spin cross-sections 
is less stringent. 
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Since the neutralino is expected to be non relativistic with average kinetic 
energy < T >w 4:0KeV{'m^/100GeV), it can be directly detected mainly via 
the recoiling of a nucleus {A,Z) in elastic scattering. In some rare instances the 
low lying excited states may also be populated [HI • In this case one may observe 
the emitted 7 rays. Finally one may be able to observe the electrons produced 
in the LSP-nucleus collision. 

In every case to extract from the data information about SUSY from the 
relevant nucleon cross section, one must know the relevant nuclear matrix ele- 
ments 50]— The static spin matrix elements used in the present work can 
be found in the literature 

Anyway since the obtained rates are very low, one would like to be able to 
exploit the modulation of the event rates due to the earth's revolution around the 
sun Q2|> assuming some velocity distribution |12I13I14( - [T5| for the LSP. One 
also would like to exploit other signatures expected to show up in directional 
experiments ||16j . Finally in a novel proposal one may be able to observe the 
reaction produced electrons |Ej instead of the standard nuclear recoils or even 
better to observe them in coincidence with the recoiling nuclei. 



2 Rates 

The differential non directional rate can be written as 

dRundir = .™ d(T{u,v)\v\ (1) 

TO^ AiriN 

where da{u,v) was given above, p(0) = 0.3GeV/cm^ is the LSP density m our 
vicinity, m is the detector mass and nix is the LSP mass 

The directional differential rate, in the direction e of the recoiling nucleus, 

is: 

dRd^r = ^-^\v\v.e 0iv.e) da(«, v 5{^^ - v.e) (2) 
AniN 27r ^rV\/2 

where 0{x) is the Heaviside function and: 

du 

where u the energy transfer Q in dimcnsionlcss units given by 



da{u, v) === 7777717X2 ii^sF{uf + Ssp^nFuiu)] (3) 



Qo 

with b is the nuclear (harmonic oscillator) size parameter. F{u) is the nuclear 
form factor and Fn (u) is the spin response function associated with the isovector 
channel. 

The scalar cross section is given by: 
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(since the heavy quarks dominate the isovector contribution is neghgible) . af^ 
is the LSP-nucleon scalar cross section. The spin Cross section is given by: 

^spin — ( 7 r) Cspin: Cspin — To ^{^^ (^) 

MHW 3(1 + 4)2 

J A 

Sin) « 5(0) = [AfioiO))' + 2^f2o{0)f2m + f2,{0)f] (7) 

J A J A 

/a I /A isoscalar and the isovector axial current couplings at the nucleon 

level [in|. 



3 Results 

To obtain the total rates one must fold with LSP velocity and integrate the 
above expressions over the energy transfer from Qmin determined by the detec- 
tor energy cutoff to Qmax determined by the maximum LSP velocity (escape 
velocity, put in by hand in the Maxwellian distribution), i.e. v^^sc — 2.84 vq with 
Vq the velocity of the sun around the center of the galaxy (229 Km/ s). 



3.1 Non directional rates 

In a previous paper jl9| we have shown that, ignoring the motion of the Earth, 
the total non directional rate is given by 



K 



Ccoh{A, ^J,r{A))ap^^o + Cspm(A, ^,.(A))cr^^_^" (spi: 



(8) 



where K 



P(0) 



y/ (v^) and 



Ccoh{A,HriA)) = 



fir (A) 



1 2 



A tcohiA) , Cspin{A, ^r{A)) — 



flriA) 



1 2 



(A) 



A 



(9) 

where t is the modification of the total rate due to the folding and nuclear struc- 
ture effects. It depends on Qmin, i-e- the energy transfer cutoff imposed by the 
detector and a = [firbvoy/2]~^ The parameters Ccoh{A, iir{A)), Cspin{A, iir{A)), 
which give the relative merit for the coherent and the spin contributions in the 
case of a nuclear target compared to those of the proton, are tabulated in table 
^for energy cutoff, Qmin = , 10 keV. Via Eq. JHJ we can extract the nucleon 
cross section from the data. 

Furthermore we have seen that ignoring the isoscalar axial current and using 
n1 = 1.22 and [2f = 2.8 for '^'^'^l and ^^F respectively we find: 



spxn 
P,X° 



{A, fir{A)) 



xlO^(A = 127), wxl02(A=19) (10) 
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Table 1. The factors cl9 = Ccoh(19, Mr(19)), sl9 = Csp»„(19, Atr(19)) and cl27 = 

Ccoh(127, ^r(127)), Sl27 = Cspin ( 127, 127)) for two values of Qmin- 



Qmin 








(GeV) 








keV 




20 30 40 


50 


60 


80 


100 


200 





cl9 


2080 2943 3589 


4083 


4471 


5037 


5428 


6360 





sl9 


5.7 8.0 9.7 


10.9 


11.9 


13.4 


14.4 


16.7 





cl27 


37294 63142 84764 101539 114295 


131580 142290 162945 





sl27 


2.2 3.7 4.9 


5.8 


6.5 


7.6 


8.4 


10.4 


10 


cl9 


636 1314 1865 


2302 


2639 


3181 


3487 


4419 


10 


sl9 


1.7 3.5 4.9 


6.0 


6.9 


8.3 


9.1 


11.4 


10 


cl27 


11660 24080 


36243 


45648 


58534 


69545 


83823 


10 


sl27 


0.6 1.3 


1.9 


2.5 


3.3 


4.0 


5.8 



It is for this reason that the hmit on the spin proton cross section extracted from 
both targets is much poorer. The form factor favors the hghter system JS] both 
for the spin and the coherent process, t(127) < i(19). In the case of the spin this 
advantage is not offset by the larger reduced mass. It is even enhanced by the 
spin ME (see Table ^1. For the coherent process, however, the light nucleus is 
no match (see Table P). 

If the effects of the motion of the Earth around the sun arc included, the total 
non directional rate is given by 

R = K [ccoh{A,fir{A))al^o{l + h{a,Q„^in)cosa)] (11) 

and an analogous one for the spin contribution, h is the modulation amplitude 
and a is the phase of the Earth, which is zero around June 2nd. The modulation 
amplitude would be an excellent signal in discriminating against background, 
but unfortunately it is very small, less than two per cent. Furthermore for in- 
termediate and heavy nuclei, it can even change sign for sufficiently heavy LSP 

CHI- 

3.2 Directional Rates. 

Since the sun is moving around the galaxy in a directional experiment, i.e. one 
in which the direction of the recoiling nucleus is observed, one expects a strong 
correlation of the event rate with the motion of the sun. In fact the directional 
rate can be written as: 

Rdir = [ccoh{A, ^r{A))a^-^o (1 + hraCOs{a - Ura Tt))] (12) 

and an analogous one for the spin contribution. The modulation now is hm, 
with a shift amTr in the phase of the Earth a, depending on the direction of 
observation. K/(27r) is the reduction factor of the unmodulated directional rate 
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Fig. 1. The expected modulation amplitude hm for A = 127 in a direction outward 
from the galaxy on the left and perpendicular to the galaxy on the right as a function 
of the polar angle measured from the sun's velocity. For angles than 7r/2 it is irrelevant 
since the event rates are tiny. 



Table 2. The parameters t, h, K,hm and am for Qmin = 0. The results shown are 
for the light systems. -\-x is radially out of the galaxy {& — 7r/2,<? = 0), +z is in the 
the sun's motion and +y vertical to the plane of the galaxy so that {{x,y,x) is right- 
handed. Om = 0,1/2,1,3/2 implies that the maximum occurs on June, September, 
December and March 2nd respectively. 



type t h dir k hm ct- 



+z 0.0068 0.227 1 
dir +(-)x 0.080 0.272 3/2(1/2) 

+ (-)y 0.080 0.210 (1) 

-z 0.395 0.060 

all 1.00 

all 0.02 



relative to the non-directional one. The parameters k , hm , ctm strongly depend 
on the direction of observation. We prefer to use the parameters k and hm, since, 
being ratios, are expected to be less dependent on the parameters of the theory. 
In the case oi A = 127 we exhibit the the angular dependence of hm for an LSP 
mass of = lOOGey in Fig.^ We also exhibit the parameters t, h, k, hm and 
am for the target A = 19 in Table |21 (for the other light systems the results are 
almost identical). 

The asymmetry is quite large in the direction of the sun's motion is large |18| . 
« 0.97. In the plane perpendicular to the sun's velocity the asymmetry equals 
the modulation. 

For a heavier nucleus the situation is a bit complicated. Now the parameters k 
and hm depend on the LSP mass as well, (see Figs El and EJ. The asymmetry and 
the shift in the phase of the Earth are similar to those of the yl = 19 system. 
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Fig. 2. The parameter k as a function of tlie LSP mass in tiie case of tlie A — 127 
system, for Qmin ~ expected in a plane perpendicular to the sun's velocity on the 
left and opposite to the sun's velocity on the right. 




Fig. 3. The modulation amplitude hm in a plane perpendicular to the sun's velocity 
on the left and opposite to the suns velocity on the right. Otherwise the notation is 
the same as in Fig|5| 



3.3 Transitions to excited states 

Incorporating the relevant kinematics and integrating the differential event rate 
dR/du from Umin to Umax wc obtain the total rate as follows: 

Rexc = ^(1 - ^)du , = ^du (13) 

where We^c = JmNQo ^^'^ excitation energy of the final nucleus, Umax = 

{y/aY-{E^/Qo) , y = v/upsilono and /Qo, Qmin (imposed by the 

detector energy cutoff) and Umax — [Vesc/aY imposed by the escape velocity 
(yesc = 2.84). 

For our purposes it is adequate to estimate the ratio of the rate to the excited 
state divided by that to the ground state (branching ratio) as a function of the 
LSP mass. This can be cast in the form: 

B RR '^ea:c(0) ^exci^Uexa '^umax^{^ ^exci^Uexc^ Umax^ COS q] fl4) 

Sgs{0) 'I^gs{Umin)[l + h{Umin) COS a] 

in an obvious notation |18p. Sgs(0) and S'ea;c(0) are the static spin matrix ele- 
ments As we have seen their ratio is essentially independent of supersymmetry, 
if the isoscalar contribution is neglected. For ^^^I it was found i9| to be be about 
2. The functions are given as follows : 

fiy/a)^ c f \ 

1^gs{Um^n)= ^^Ff^^) [V'(a^/^) - V^e.c)] (15) 
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[ipia^/uil + Uexc/u)) -^^:{yesc)] du (16) 

The functions -0 arise from the convolution with LSP velocity distribution. The 
obtained results are shown in Fig. ^ 




rriLSP 



[GeV) 




rriLSP 



(GeV) 



Fig. 4. The ratio of the rate to the excited state divided by that of the ground state as 
a function of the LSP mass (in GeV) for ^^^/. We assumed that the static spin matrix 
element of the transition from the ground to the excited state is a factor of 1.9 larger 
than that involving the ground state, but the functions Fii(u) are the same. On the 
left we show the results for Qmin = and on the right for Qmin = 10 KeV. 



3.4 Detecting recoiling electrons following the LSP-nucleus 
collision. 

During the LSP-nucleus collision we can have the ionization of the atom. One 
thus may exploit this signature of the reaction and try to detect these electrons 
[T7) . These electrons are expected to be of low energy and one thus may have a 
better chance of observing them in a gaseous TPC detector. In order to avoid 
complications arising from atomic physics we have chosen as a target ^"^Ne. Fur- 
thermore, to avoid complications regarding the allowed SUSY parameter space, 
we will present our results normalized to the standard neutralino nucleus cross 
section. The thus obtained branching ratios are independent of all parameters 
of supersymmetry except the neutralino mass. The numerical results given here 
apply in the case of the coherent mode. If, however, we limit ourselves to the 
ratios of the relevant cross sections, we do not expect substantial changes in the 
case of the spin induced process. 

The ratio of the our differential (with respect of the electron energy) cross 
section divided by the total cross section of the standard neutralino-nucleus 
elastic scattering, nuclear recoil experiments (nrec), takes .17 the form: 



'■\F(uM£ + AW 

,kdT. (17) 



S\ da £ dS,K^-^[F{pM^ + A))f 
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where 



K = , K = ^ , a = P^.k, e = q.K, 

Px Px 



— Y + — 1)] and 2^p^^ = 2/^^^^ is the momentum q trans- 

ferred to the nucleus and F{q) is the nuclear form factor. The outgoing electron 
energy lies in the range < T < -^T^ — Sni- Since the momentum of the out- 
going electron is much smaller than the momentum of the oncoming neutralino, 
i.e. K K 1, the integration over ^ can be trivially performed. 

We remind the reader that the LSP- nucleus cross-section anrec takes the 
form: 



^J^fA'ap C 2dm'^^^rvOY (18) 

^^r[JP) Jo 

ase of ^°Ne he bindi 

by [HI: 



'Hr{p)' 

In the case of ^°Ne he binding energies and the occupation probabilities are given 



e„i = (-0.870, -0.048 , -0.021) , = (2/10, 2/10, 6/10). (19) 

in the obvious order: (ls,2s,2p). In Fig [S] we show the differential rate of our 
process, divided by the total nuclear recoil event rate, for each orbit as well as 
the total rate in our process divided by that of the standard rate as a function of 
the electron threshold energy with threshold energy in the standard process. 
We obtained our results using appropriate form factor [TT] . 

From these plots we see that, even though the differential rate peaks at low 
energies, there remains substantial strength above the electron energy of 0.2 keV^ 
which is the threshold energy for detecting electrons in a Micromegas detector, 
like the one recently (20] proposed. 



T 



T keV 




25 1.5 1.75 E, 



Fig. 5. Shown on the left is the differential rate, divided by the total rate associated 
with the nuclear recoils, as a function of the electron energy T (in keV). Each atomic 
orbit involved in the target ^°Ne is included separately. The full line, the short-dashed 
line and the long-dashed line correspond to the orbits Is , 2s and 2p respectively. 
Shown on the right is the ratio of the total rate for the novel process divided by that 
of the standard process as a function of the electron threshold energy, assuming zero 
threshold energy for the standard process. This ratio may increase if such a threshold 
is included. 
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4 Conclusions 

Since the expected event rates for direct neutralino detection are very low |5l8) . in 
the present work we looked for characteristic experimental signatures for back- 
ground reduction, such as: 

• Standard recoil experiments 

Here the relevant parameters are t and h. For light targets they are essentially 
independent of the LSP mass 2H1> essentially the same for both the coherent 
and the spin modes. The modulation is small, h « 0.2%, but it may increase 
as Qmin increases. Unfortunately, for heavy targets even the sign of h is 
uncertain for Qmin = 0. The situation improves as Qmin increases, but at 
the expense of the number of counts. 

• Directional experiments ^B] Here we find a correlation of the rates with the 
velocity of the sun as well as that of the Earth. One encounters reduction 
factors k/2tt, which depend on the angle of observation. The most favorable 
factor is small, w l/47r and occurs when the nucleus is recoiling opposite 
to the direction of motion of the sun. As a bonus one gets modulation, 
which is three times larger, hm =~ 0.06. In a plane perpendicular to the 
sun's direction of motion the reduction factor is close to l/127r, but now the 
modulation can be quite high, hm ~ 0.3, and exhibits very interesting time 
dependent pattern (see Table 13 Further interesting features may appear in 
the case of non standard velocity distributions |15) . 

• Transitions to Excited states 

We find that branching ratios for transitions to the first excited state of 
^^^I is relatively high, about 10%. The modulation in this case is much 
larger he^c ~ 0.6. We hope that such a branching ratio will encourage future 
experiments to search for characteristic 7 rays rather than recoils. 

• Detection of ionization electrons produced in the LSP collision 

Our results indicate that one can be optimistic about using the emitted 
electrons in the neutralino nucleus collisions for the direct detection of the 
LSP. This novel process may be exploited by the planned TPC low energy 
electron detectors. By achieving low energy thresholds of about 0.25 keV, 
the branching ratios are approximately 10 percent. They can be even larger, 
if one includes low energy cutoffs imposed by the detectors in the standard 
experiments, not included in the above estimate. 

As we have seen the background problems associated with the proposed 
mechanism are not worse than those entering the standard experiments. In 
any case coincidence experiments with x-rays, produced following the de- 
excitation of the residual atom, may help reduce the background events to 
extremely low levels. 
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